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INTRODUCTION 


The use of linear elastic fracture mechanics In the design of structures 
requires the knowledge of the fracture toughness of the material. For this 
purpose the ASTM has established a standard test method. (E399) for obtaining 
the plane-strain fracture toughness of a given material, [1]. One of the 
standard specimens used in such a test is the compact tension specimen (CTS) 
shown in figure 1. 

Very rigid standards have been Imposed in reference [1] to insure that 
valid fracture toughness values will be obtained in a given test. One of the 
requirements is the production of a fatigue crack. However, the fatigue crack 
will usually not grow uniformly across t e specimen thickness, i.e. the initially 
straigth crack front will become curved. The crack length thus varies across 
the specimen thickness and some average value must be used to calculate the 
fracture toughness. A standard measure is provided in reference [1] such that 
if the crack front curvature, as measured after the specimen is broken, is great- 
er than this standard measure, the test is invalid. 

Recently, however, suggestions have been made, [2] that the above standard 
is too rigid and can be relaxed since the effect of the curved crack front may 
be less than originally anticipated. These suggestions, however, are based 
on rather weak evidence, [3] and it is not yet clear whether the E 399 standard 
should be relaxed or not in this respect. What Is needed is a clear definitive 
analysis of the effect of the crack front curvature on the stress intensity 
factor. The objective of this Investigation is to perform this analysis. 

Several three-dimensional finite element analyses have been made of the 
compact tension specimen, e.g. [4,5]. These analyses showed that for the CTS, 
the stress intensity factor is a maximum at the center of the specimen and a 
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minimum at the surface of the specimen. When crack curvature is introduced, 
the maximum stress intensity factor moves to the surface of the specimen and the 
minimum occurs at the center. The dimensions of the specimens treated in these 
analyses differed from those of the standard specimen, A few attempts have also 
been made to determine experimentally the effect of crack tip curvature using 
a 3-D stress freezing technique [6], and scattered light speckle interferometry 
[7]. Both these Investigations showed qualitative agreement with the above 
results. 

The present investigation uses the method of lines (MOL) to obtain an 
accu'ate three-dimensional solution for the standard CTS with "urved crack 
fronts. The MOL is described fully in references [8] and [9], where it is also 
shown that accurate solutions of crack problems can be obtained. An improved 
version of the method and the application to the standard CTS with a straight 
crack front has been presented in a previous paper, [10]. The method will be 
only briefly reviewed here for completeness. 

BACKGROUND OF THE METHOD OF LINES 

Within the framework of linearized elasticity theory, the equations of 
elastic equilibrium in terms of displacements are 

G vV + (X + G) 0 * 0 , i = 1, 2, 3 (1) 

1 >1 

9 “ u^ . 
j.3 

where X and G are Lame’s constants, u^ are the displacement and the usual tensor 
notation is used. 

For a finite geometry solid we construct three sets of parallel lines, 

Fig. 2. Each set of lines is parallel to one of the coordinate axes and thus 
perpendicular to the corresponding coordinate plane. An approximate solution 
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of Eq. (1) can then be obtained by developing solutions of ordinary differential 

equations along tlie x-directional lines. As seen in the figure, there are a 

total of Jl ■ NYxNZ such lines where NY is the number of lines along the 

y-direction and NZ is the number of lines along the z-direction in a given 

plane, respectively. We define the displacements along these lines as u^^, u^, 

. . . u^. The derivatives of the y-directional displacements on these lines 

with respect to y are defined as v' [ . v' | 8"^ the derivatives 

12 ^ 

of the z-directional displacements with respect to z are defined as w' I » w' I , 

1 2 

. . ., w'l^. These displacements and derivatives can then be regarded as 
functions of x only since they are variables on x-directional lines. When these 
definitions are used, the ordinary differential equation along a generic line 
ij using a 5-point difference formula, may be written as: 


^ . g - 2v) , 

,2 2(1 - V) 4zh 2 ^ i-2,j 

dx y 




30u^j - 16u 


i+l,j 


+ J + 




*i+2,j" ' 12h;.2 '“i,j-2 ‘ ^^^i,j-l “ ^^^i,j+l 




i,j+2"-' dx 2(1 - v) 


( 2 ) 


where 


f. .(x) * v' + w'; v' = ^ and w’ « — 
ij dy dz 


(3) 


Sinilar differential equations are obtained along the other x-directional lines. 

Since each equation has the terms of the displacement on the surrounding lines. 

These equations constitute a system of ordinary differential equations for the 

displacements u , u . . . , u„ . It is to be noted that for lines lying on or 
12 *- 

near a boundary appropriate noncentral difference formulas must be used. Also, 
the shear boundary conditions are incorporated in these equations as described 
in reference [9]. 
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The set of i. second order differential equations represented by Eq. (2) 

* 

can be reduced to a set of 2i first order differential equations by treating 
the derivatives of the u's as an additional set of I unknowns. 


u 


U1 




dT • 


W 


The resulting 21 equation*: can now be written as a single first order matrix 
differential equation 


^ ■ A U + 
dx X 


(5) 


where 1) and R are column matrices of 22. elements each and Aj is a 22 x 22 matrix 
of the constant coefficients appearing in Eqs. (2) and (4). 

In a similar manner, to solve Eq. (1) for 1=2 and 3 ordinary differen- 
tial equations are constructed along the y- and z-directicnal lines, respectively. 
These equations are also expressed in an analogous form to Eq. (5); they are 


dy 

dz 


A V + 
2 

A W + 

3 


dS(y) 

dy 

dT(z) 

dz 


( 6 ) 

(7) 


Equations (5) to (7) are linear first-order ordinary matrix differential 
equations. They are, however, not independent, but are coupled through the 
vectors R, S, and T whose components are given by equations similar to Eq. (3). 
They must therefore be solved iteratively correcting R, S and T after each 
iteration. The elements of the coefficient matrices A , A , and A are all 
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constants, being functions of the coordinate increments and Poisson's ratio 
only. 


To solve equations (5) - (7) a recurrence relation algorithm was set up 
as described in detail in reference (10). This algorithn' is particularly useful 
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for solving two-point boundary value problems, requiring one sv/eep through from 
one boundary to the opposite one to obtain the missing boundary conditions and 
a second sweep through to obtain the complete solution everywhere. 

CALCULATION OF STRESS INTENSITY FACTOR » 

The stress intensity factor is defined as follows: 

Kj = liim a ( 2 ttR)" 

R->0 ^ 


where R is measured from the crack tip and is normal to the crack front, n Js 
the singularity. It was found, however, that due to the coarseness of the grid 
used, the usual plotting and extrapolating techniques gave results that were 
erratic and of questionable accuracy. This was compounded by the fact that the 
precise crack tip location is not really known except that it is approximately 
midway between two lines, one of which has zero stress specified. It was found, 
however, that by using two terms in the stress and displacement series expansions 
around the crack tip, good revolts could be obtained even with the coarse grid 
used. Furthermore, this also permitted the determination of the actual crack 
tip location from the computed results. The method utilized is as follows, [11], 
we take 


''iy-O ■ 




+ r 

2tt 


L, 

+ V(R + r)’ 


( 8 ) 


'^y [y«0 


K. [ ^ r 

/2tt(R - r) I 


(9) 


where a is a function of Poisson's ratio, n was assumed to be -hy and r is the 
crack edge position corr ction measured from the originally assumed midpoint 
position. Using displacement data from three adjacent nodes to the crack edge 


5 



OF POOR QUAUITT 


in Eq. (8), values ofaK L /K_ and r are calculated for ec-ch value of z, 

i , 1 i 

with R also measured frcm the halfway point between nodes specifying boundarv 
stresses and displacements, respectively. Note that a would be equal to 3.56 
for the plane strain case and A.O for the plane stress case. In the case of 
curved crack fronts a normal is drawn at the point where Is determined. Tlie 
displacement values at the points where the normal meets the grid lines are extra- 
polated using the nodal displacement values. These extrapolated values of 
placement are used to obtain the local K^. values for the curved ere :k froi 

RESULTS AND DISCUSSION 

The MOL was first applied to a compact tension specimen with tensile load- 
ing. The specimen dimensions were W/a = 2.0, B/a =• 1.0 and L/a ■ 2.4. Wiiere 
a, B, W and L are the crack length, specimen thickness, w. dth and length respec- 
tively. Crack length to width ratio was maintained close to 0.5. Figures 2 
and 3 show the coordinate system and line number designation. 

Figure 4 shows a plot of the local stress intensity factor variation 
through the thickness of the specimen with tensile loading. Curve 1 is for a 
straight crack front. For this case the maximum SIF occurs at the center of 
the specimen. It drops by 13 percent on the surface. This trend is in accor- 
dance with the previously reported results of Raju (12). 

The remaining curves in the figure are for increasing crack tunnel depth. 
Crack tunnel depth is defined as the difference between the crack lengths at 
the center and the surface of the specimen. Since approximate parabolic curves 
were used in the study, it was difficult to assign a single crack curvature for 
the crack front. An increasing crack tunnel depth reflects an increasing crack 
curvature. For crack tunnel depth of 0.1 we see that SIF decreases at tiie center 
by 9.5 percent, while at the surface its value increases by 20.5 percent. Simi- 
lar type of variation is reported by Pereira [4]. However, a direct comparison oi 
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the results Is not possible^ because a different crack and specimen size were 
used in that study. Herein, the dimensions prescribed in ASTM rtandard [1] 
are used. 

For higher crack tunnel depths tbe maximum SIF value no longer occurs at 
the surface, but shifts to the adjacent interior point of the grid. This 
trend is not reported in reference [4] . There could be two possible reasons 
for this differece, 1) the crack length a/W used in reference [4] is 0.25 
which is smaller than the one used in the present study, 2) three layers of 
elements were taken in the thickness direction, possibly causing the maximum 
SIF to be missed. 

Figure 5 shows the \ariation of the non-dimensional SIF for the same 
specimen with a parabolic shear load along the load line. The results are 
similar to those of fig'i.e 4 for the tensile loading. 

Figure 6 shows the variation of center (r), surface (T.) and thickness 
average (av) stress intensity factors with increasing crack tunnel depth for a 
standard compact tension specimen under shear loading. The stress intensity 
factors are non-dimensionalized with the respective values of SIF occur- 

ing for a straight crack front (SC). We find that non-dimensional SIF value at 
the center of the specimen keeps on decreasing as the crack tunnel depth in-- 
creases. However, the value of SIF at the surface increases to a maximum and 
tlien decreases. The thickness average SIF constantly decreases with increasing 
crack tunnel depth. Similar trend for the thickness average SIF was observed 
by McGowan [6] in his study on single edge notch specimen. This figure shows 
clearly that the presence of a curved crack front can significantly change the 
value of the stress intensity factor. 

In figure 7 the change in the thickness average SIF is plotted with a non- 
dimensional factor n. The value of n corresponds to the difference between the 
average crack length and the surface trace of the crack divided by tiie 
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average crack length (a ) . The calculation of a la done by coaputlng an 
average of the three c.ack measurenents taken at the center of the crack front 
and midway between the center and the end of the crack front on each side, which 
Is the definition used in [1]. Reference [1] restricts the length of either 
surface trace of the crack to greater than 90 percent of the average crack length 
(a }. Based on this criterion, the ASTM region la narked on the figure 7. Any 

o V 

test in which the specimen has a surface trace within the marked region is con- 
sidered a valid test. For the worst case when the surface trace is 90 percent 
of a in the specimen, there is a 7 percent lowering in the thickness average 
stress intensity factor, which indicates that the use of plane strain formula 
given in reference [1} will result in an overestimation of the fracture toughness 
of a material by 7 percent. 

On the basis of the experimental results and three-dimensional analysis of 
Neale and Pereira [4,5], an amendment has recently been Incorporated, [13] to 
extend the limit on the surface trace to 85 percent of average crack length. 

This figure corresponds to a n value of 15 percent. From the figure 7 we find 
a lowering of 12.5 percent in the thickness average SIF for this value of n* 

This implies that the use of two-dimensional plane strain formula will result 
in an overestimation of the fracture toughness of a material by 12.5 percent. 
These figures thus show the effect of crack curvature on the overestimation of 
the fracture toughness for the compact tension specimen. 

CONCLUSIONS /iND RECOMMENDATIONS 

1) The results from the niunerlcal analysis show that the stress intensity factor 
decreases at the center of the compact specimen due to crack curvature. For 
the low values of crack tunnel depth the maximum stress intensity factor is 
found to be at the surface of the specimen. If the crack tunnel depth is 
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Increased, the tnaxlmtan value of stress intensity factor no longer remains 
at the surface, but shifts to an interior point of the specimen, stil] 
remaining close to the surface. 

2) The average value of stress intensity factor through the thickness keeps on 
decreasing with increasing crack tunnel depth. 

3) The ASTM test method LlJ allows either surface trace of the crack to be less 
Chan 90 percent of the average crack length. Based on an average SIF calcu- 
lation, this will overestimate the fracture toughness by 7 percent. It is 
further estimated that the new amendment changing the 90 percent to 85 per- 
cent will lead to an overestimation of fracture toughness by 12.5 percent. 

The piots for stress Intensity factor through the thickness of compact 

tension specimen for two different types of loading predict that the crack growth 
should start at the surface of the specimen first. This conclusion from analyti- 
cal study is contrary to the experimentally observed crack propagation. It is 
anticipated that plastic flow can significantly change the stress conditions 
existing at the surface, which is, of course, not accounted for in an elastic 
analysis. Neale [14] developed an elasto-plastlc model to fix this discrepancy, 
but it is an approximate analysis. To fully study the problem, we need to 
analyze a full three-dimensional elastic-plastic case. The results of such a 
study will be presented in a subsequent publication. 
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Figure 3. Coordinate system and line numbering for 
single edge notch and compact tension 
specimens. 
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Figure.4. Vsiriation of dimensionless SIF with thickness for different 
c' ack tunnel depths for compict tension specimen under 
tensile loeding. 
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Fifure 5. Variation of dimensionless SIF through the thickness for 
different crack tunnel depths, for compact tension 
specimen under paraholically applied shear load. 
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Figure 6. Variation of canter, surface and thickness 
average stress intensity factor with 
increasing crack tunnel depth. 
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Figurs 7. Percent variation in average stress intensity 
factors between straight crack and curved 
crack fronts. 
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